Bound, virtual and resonance S'-matrix poles from the Schrodinger equation 
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A general method, which we call the potential 5-matrix pole method, is developed for obtaining 
the 5-matrix pole parameters for bound, virtual and resonant states based on numerical solutions of 
the Schrodinger equation. This method is well-known for bound states. In this work we generalize 
it for resonant and virtual states, although the corresponding solutions increase exponentially when 
r — > co. Concrete calculations are performed for the 1 + ground and the + first excited states of 
14 N, the resonance 15 F states (l/2 + , 5/2 + ), low-lying states of 11 Be and n N, and the subthreshold 
resonances in the proton-proton system. We also demonstrate that in the case the broad resonances 
their energy and width can be found from the fitting of the experimental phase shifts using the 
analytical expression for the elastic scattering .S'-matrix. We compare the S'-matrix pole and the 
_R-matrix for broad swa resonance in 15 F. 

PACS numbers: 26.20.+I, 24.50.+g, 25.70.Ef, 25.70.Hi 



I. INTRODUCTION 

Analysis of the S'-matrix pole structure is a powerful 
method in quantum physics. It is well-known that the 
poles of the S'-matrix in the complex momentum (or en- 
ergy) plane correspond to bound, virtual and resonance 
states. There is a well-known relation between the S- 
matrix and Jost functions, the singular solutions of the 
Schrodinger equation at r — ► 0. The conventional nu- 
merical method for bound states is to search for solu- 
tions, which only have an outgoing wave at pure imagi- 
nary momenta in the upper half momentum plane. The 
corresponding wave function is an exponentially decreas- 
ing solution when r — ► oo. 

Virtual or resonance states are described by the wave 
functions containing only the outgoing waves asymptot- 
ically, which exponentially increase due to the complex 
momenta. In the past (see [l[ and references therein) , the 
analytical continuation onto the unphysical energy sheet 
of the Lippmann-Schwinger as well as the momentum- 
space Faddeev integral equations were used to find the 
resonance properties. The normalization formula for the 
bound state vertex function in the momentum space was 
generalized in Q for the resonance and virtual states. 

Such states are considered as unaccomplished bound 
states. This means that a bound state pole should move 
down the positive semi-axis of the complex momentum 
plane while the interaction strength decreases. At some 
critical value of the interaction strength, the pole ap- 
proaches the zero energy point, which belongs to the 
contour of integration. After a subsequent decrease of 
the interaction strength, the pole moves to the lower 
part of the complex momentum plane (unphysical energy 
sheet) dragging the integration contour in the Lippmann- 
Schwinger equation to secure the convergence of the in- 
tegral. This leads to the appearance of an extra term in 
the right-hand-side of the equation, which is the residue 
of the integrand at the pole. 



This method of the analytical continuation has been 
applied successfully to different physical systems. Unfor- 
tunately, it can not be used directly in the case of charged 
particles. We should also note that an analytical form of 
the Fourier transform of the potential, which is an input 
in the Lippmann-Schwinger integral equations, is known 
only for a limited number of potentials. 

The problem of the exponential increase of the Gamow 
resonance wave function in the asymptotic region can be 
solved by a complex scaling method based on the so- 
called ABC— theorem [3J . This method consists of solving 
the Schrodinger equation on a ray in the first quadrant of 
the radial complex plane rather than on the real axis of 
the coordinate r. This ray can be obtained by the follow- 
ing transformation of the radial coordinate r and the con- 
jugate momentum p: r — > r exp(i9) and p — > p exp(— i#). 
As a result, the bound state spectrum is supplemented by 
the S-matrix poles situated in the sector defined by the 
angle 9 between the real axis and the ray in the fourth 
quadrant of the complex momentum plane. The axis ro- 
tation angle, 9, is limited by the position of the potential 
singularities in the radial complex plane. It is important 
that the complex scaling method can be applied to the 
case of charged particles. The method is valid because 
the Coulomb potential satisfies the scaling condition of 
the ABC— theorem. An application of this method to 
resonances in nuclear reactions was presented in Q ■ The 
numerical realization of this method is a rather complex 
one. 

A few different techniques to determine the resonance 
energy, width and resonance wave function based on the 
solution of the Schrodinger equation have been previously 
suggested. In these methods the normalization of the 
resonant wave function is achieved using the Zel'dovich's 
normalization [Bj], which is difficult in practical realiza- 
tion due to slow convergence of the integrals. First we 
refer to the method of solution of the radial Schroodinger 
equation to determine resonances suggested in [6| . In this 
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method the complex eigenvalue and the Gamow wave 
function can be found by integration of the Schrodinger 
equation imposing the boundary conditions in the ori- 
gin and the asymptotic region. To solve the equation 
the Fox-Goodwin numerical method was applied and the 
logarithmic derivatives of the internal and external wave 
functions were matched. However, this method fails in 
the vicinity of the threshold, for broad and subthreshold 
broad resonances (imaginary part of the momentum is 
larger than its real part) and antibound (virtual) states. 
We underscore also that this method can be applied only 
for the potentials with the finite interaction radius be- 
cause of the problem with numerical calculation of the 
exponentially increasing wave function. The application 
of the method for the unstable nuclei can be found in 
@,|. 

A pole search has also been used in 9] by solu- 
tion of the Schrodinger equation with the short range 
interaction for the scattering wave function. The 
Zel'dovich's normalization procedure for the Gamow res- 
onance wave function supplemented by the exterior com- 
plex scaling [10] was used. The norm of the Gamow reso- 
nant wave functions does exists for charged particles also 
(ToL [Tl| . The method allows one to find resonances and 
even subthreshold resonances but it cannot be applied to 
the virtual states. 

The method^ which is also close to our approach, was 
discussed in The asymptotic wave function in this 
method contains auxiliary S-matrix which coincides with 
the physical S-matrix at the resonance poles at which the 
solution becomes pure outgoing wave. The method was 
applied for determination of the low-energy 5 He and 5 Li 
resonance parameters (12J. 

In the present work, we demonstrate how to find the 
poles of the S-matrix corresponding to bound, virtual 
and resonance states and the S'-matrix residues in these 
poles by solving the Schrodinger equation with the nu- 
clear plus Coulomb potentials using the analytical prop- 
erties of the S'-matrix. In contrast to the previously pub- 
lished methods, in our S-matrix pole method the nor- 
malization of the resonant wave function is based on 
the connection between the residue of the S-matrix in 
the pole and the asymptotic normalization coefficient 
(ANC). This relationship is universal and can be ap- 
plied to bound, virtual and resonance (narrow and broad) 
states (Tl| making our technique universal, and that is 
the main distinction of our method from the previously 
published ones. The ANC is the amplitude of the tail 
of the bound, virtual or resonant wave function [TU . [l3| . 
For the resonant state, the ANC is related to the reso- 
nance width [ll| . The use of the ANC doesn't require the 
normalization of the state corresponding to the S-matrix 
pole and this why our method allows one to determine 
both narrow and broad resonances, and even antibound 
states. 

A simple relation between the ANC (nuclear vertex 
constant) and the overall normalization of the peripheral 
astrophysical S-factor suggested in [lj, [l5|, [lj| makes it 



extremely important for obtaining astrophysical S fac- 
tors. Note that the normalization method proposed by 
Zel'dovich [j| was generalized in [ll[ for the interaction 
potential with a Coulomb tail. 

The S-matrix pole method addressed here has been ap- 
plied earlier to the virtual (singlet) deuteron and virtual 
triton with different short-range potentials. The results 
of the two-step Gamov state normalization for the vir- 
tual (antibound) state of 3 H were published in [T3|. For 
charged particles, the virtual state becomes a subthresh- 
old resonance [18|]. Here we present new results for the 
subthreshold resonances for the ground state of 2 H. We 
also calculate the ground and the first excited states of 
14 N and the resonance states of 15 F. Finally, our method 
is applied to the three lowest T = | states in 11 Be and 

n N. Considering the ^ + state in n N as an example, we 
demonstrate how to determine in a model-independent 
way the energy and width of the broad resonance us- 
ing the S-matrix analytical structure, which includes the 
resonant pole. Moreover, we demonstrate that the po- 
tential S-matrix pole method addressed here gives the 
resonance energy and width, which are very close to the 
model-independent results obtained from the analytical 
expression for the S-matrix in the vicinity of a single pole 



We use the system of units in which h=c=l. 

II. A NUMERICAL CALCULATION METHOD 

To describe the nuclear interaction we adopt the 
Woods-Saxon potential 
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V N (r) = - [V Q - V LS (L ■ S)-^ — 



r> (!) 



'ml rdr 1 + exp [l^^L] 

where Vo (Vls) is the depth of the central (spin-orbital) 
potential; L is the orbital momentum operator for the rel- 
ative motion of the particles; S is the spin operator; 
is the pion mass; a is the diffuseness and Rn = tqA 1 ^ 3 
(ro is the radius parameter of the nuclear potential, A 
is the atomic mass number). The Coulomb interaction 
potential is taken in the form 

( Z 1 Z 2 e 2 ^-g _ 



Vc{r) = 



2R C 
Z 1 Z 2 e 2 



r < R c , 



(2) 



r > R c , 

where Z±e and Z^e are the charges of the particles; Rc — 
rc A 1 / 3 (rc is the parameter of the Coulomb radius). 

The radial wave function ui(r) for the partial wave 
with the orbital momentum / is the solution of the radial 
Schrodinger equation (\i\2 is the reduced mass, E is the 
energy in CM system) 

-J + 2 Ml2 [E - V (r)] - U (r) = 0. (3) 

ar z r z 



Here, ui(r) satisfies the standard boundary condition at 
the origin: 



Mr)\ r =o = °- 



(4) 
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To write the boundary condition for the derivative of 
ui(r), we analyze the behavior of the wave function 
near the origin. The sum of the potentials V(r) = 
Vn(t) + Vc(r) multiplied by r is limited. Therefore 
we choose the point tq near the origin, and denote 
k 2 = 2p, 12 [E-V(r )}. 

The solution of the Schrodinger equation 



\dr 2 



1(1 + 1) 



ui{r) = 0, 



which satisfies the condition (|U), is proportional to the 
function gi(kor) = kgr-j^kyr), where ji(kor) is the spher- 
ical Bessel function. Taking this into account, one can 
use the initial condition for Eq. © as follows 



u i( r )\r=r = 9i{kor ), u'i(r)\ r = ro = k gi(k ro). (6) 

Note that the energy E is negative for bound and virtual 
states and complex for resonance states. In the external 
region r > Rq, where the nuclear potential can be omit- 
ted with reasonable accuracy, the general solution of Eq. 
© is given by 



ui(r) 



s (r) = Cl~\k)u\~\kr) 



C, 



(+) (fc)uj +) (£x), 



(7) 

where k = y^S/Iiii?, are the coefficients that can be 
found by matching ui(r) to the solution in the internal 

region at r = R [55]. The functions ui (p) can be 
written in terms of the regular Fi(rj,p) and the irregular 
Gi(r],p) Coulomb wave functions 



ul ± \p) = e* s ? 



[G l (r ] ,p)±iF l (r ] ,p)}, 



(8) 



where r\ = Z^Z^e 2 pvilk is the Sommerfeld parameter, 8f 
is the Coulomb phase shift given by <5p = argr(l + Z + irf) 

and p = kr. The asymptotic forms of ui (p) are given 
by 



sure that c\ \k)/Ci + \k) goes to zero, we calculate the 
ratio of the Schrodinger equation solution (for the sum 
of a nuclear and the Coulomb potentials) and the out- 
going wave in the Coulomb potential. This ratio must 
approach a constant in the asymptotic region. We also 
check the equality of the logarithmic derivatives of ui (r) 

and u\ + \kr) at r — Rq- Rq should be chosen a little 
larger than the radius of the nuclear potential. Accord- 



(5) ing to the scattering theory [2l|, the vanishing of C ( 



(-) 



at 

the positive imaginary semi-axis in the complex momen- 
tum plane corresponds to the bound state, while that on 
the negative imaginary semi-axis corresponds to the vir- 
tual (antibound) state. The resonant state is defined by 
the zero of C\ ^ in the fourth quadrant of the complex 
momentum plane. 

The S'-matrix is the ratio C, (+) (fe) / 'Cf -) (fc) , which has 

a pole at k = k if Cj" ] (k ) = [20]. For the poles of S'- 
matrix of the first order the residue at the pole ko should 
be 



Res (S l (k )) = A l (k ) 



q (+) (fco) 

cl-\k y 



(ii) 



where C, (ko) is the derivative at the pole k — ko. To 

find Ai (fc ) , we calculate c\~^ (k) close enough to the pole 
ko. Then, we use the fit function 



C ( -\k) = ai (k - k ) + a 2 (k - k ) 2 



(12) 



to obtain the coefficients of the expansion a± and a 2 for 
which Ci ^'(ko) = 01. The described method we call the 
potential ^-matrix pole method. 

III. RESULTS 
A. The bound states of 14 N 



,(+) 



(p) -> exp 



p - T) In 2p - — 



00, 



(9) 



exp 



p — -q In 2p 



In 
~2 



P - 



00. 



The coefficients C ; (+) (fc) and C t 



(10) 

are proportional 



to the corresponding Jost functions [20|, |21[ . The func- 
tions flH]) and (fTU)) describe outgoing and incoming waves, 
respectively. We can solve the Schrodinger equation nu- 
merically and search for the energy at which the coef- 
ficient Cj~\k) vanishes. This condition (Cj (fe) = 0) 
means that we are dealing with only the outgoing wave 
in the asymptotic region (r — > 00). Note that for vir- 
tual and resonance states the first term in Eq. {7} is 
much smaller than the second one, which makes it dif- 
ficult to obtain a solution and an eigenvalue. To make 



To show how the method works, we start from its ap- 
plication to the bound states of 14 N considering it as a 
two-body bound state 14 N = 13 C + p. We assume that 
the proton in the lpi/2 orbital is coupled to the l/2~ 
ground state of 13 C to form the 1 + ground state and + 
excited state of 14 N. 

To describe these states in the two-body (core + mi- 
cleon) approach, we choose the geometrical parameters of 
the bound state Woods-Saxon potential to be ro = rc = 
1.2 fm and a = 0.5 feu. The well-depth procedure provid- 
ing the experimental binding energy leads to Vo = 51.65 
MeV and V LS = 1.5 MeV for the 1+ state and V = 47.71 
MeV and V LS = 1.5 MeV for the 0+ state. The coeffi- 
cients c\ + \k) and c[ \k) are found from the set of 
equations Ui{r\) — w° s (ri), u/(r 2 ) = u" s (r 2 ) (uf s is the 
known asymptotic solution) , where both the neighboring 
points 7*1 and r 2 should be chosen in the asymptotic re- 
gion. In this work we choose as an example r± = 0.5R max 
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and r 2 = OMR max , where R ma x = N R N . The pa- 
rameter N should be big enough to fulfill the condition 
ui(ri)/uf s (ri)= const. In Figs. [I]and[2j the wave func- 
tion for the 1 + state of 14 N and the ratio of the wave 
function to the Whittaker function describing its asymp- 
totic behavior are shown. For the + state, the wave 
function and its ratio to the Whittaker function are very 
similar to those of the 1 + state. From these figures, one 
can conclude that the coefficient Cj (fc) is equal to zero 
and the wave function coincides with its asymptotic form 
when r > Rq. 

From Eq. ([7]) we get that in the external region the 
radial bound state wave function is given by 

„(»-)«»( r ) =C r(+)( fe ) u (+)( jfer ) ) (13) 

where k = iKb s (k^ > 0). Normalizing the bound state 
wave function to unity we can rewrite its asymptotic term 
as 

u\ bs)as (r) = bi W_ w+1/a (2K fc8 r) 1 (14) 

where 6/ is the single-particle ANC, W_ Vbs ^ +1 / 2 {2 Kf, s r) 
is the Whittaker function determining the radial shape 
of the bound state wave function, rjb s is the Coulomb 
parameter for the bound state, Kb s = y/2 ]Ivz£bs is the 
bound state wave number and Sb s is the binding energy 
of the bound state. For the adopted geometrical param- 

— 1/9 

eters, we get &i( 9r ) = 4.250 fm ' for the ground state 

i/2 

and b\t exc \ — 2.457 fm ' for the excited state. Note 
that the single-particle ANC is sensitive to these param- 
eters [Hj]. 

The residue at the bound state pole of the S-matrix is 
given by 

Ai(k ) = (-l) l+1 ibf. (15) 



Our calculated residues of the S'-matrix at poles re- 
lated to the ground and excited states give the values 
A gr — 18.061 i fm -1 and A exc = 6.039 i fm -1 , respec- 
tively. Found from these residues, the single-particle 
ANCs coincide with &i( sr ) and &i( eX c) given above and 
found from the bound state wave functions. This vali- 
dates method of calculation of the residue of the S'-matrix 
at the bound state pole presented here. 

B. Virtual (antibound) state 

Here we apply our method to obtain the energy of the 
virtual (antibound) state in the np system at I = 0, 
taking into account only the short-range Yukawa nuclear 
potential Vjv(r) = Vbr -1 exp(— r/R). The virtual state 
corresponds to k — — i n (n > 0), i.e. the pole of the 
S-matrix is located on the negative imaginary semi-axis 
in the complex momentum plane. It generates the expo- 
nentially increasing term Uq (Ax) when r — > oo while the 

second term Uq (kr) becomes exponentially small mak- 
ing it very difficult to determine the energy (momentum) 

when Cq ^ = 0, which is the condition for the virtual 

pole. For this reason, to calculate Cq~\— in), one should 
obtain a solution with very high precision. In our calcu- 
lations, the energy | e v (N) | of the virtual state calcu- 
lated as function of R max — N Rm decreases smoothly 
as N increases. However, when N > N max the energy 
exhibits a sudden change to the larger value. It means 
that for r > R m ax the solution is not precise enough to 
calculate Cg '(—in) accurately. That is why we adopt 
(np = £v(N max ) as the virtual pole energy. Our result 
£«p = —0.067 MeV agrees very well with the one obtained 
using the integral equation method [13, H3| • The calcu- 
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FIG. 1: Comparison of the normalized radial bound state FIG. 2: Ratio of the calculated radial bound state wave 
wave function for the 14 N (1 + ) state (solid line) with the cor- function to the Whittaker function (W_ v j +1 / 2 {2K,r)) for the 
responding asymptotic form (W_ J) !+1 / 2 (2Kr), dashed line). 14 N (1 + ) state. 
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lated residue of S'-matrix in pole is A n 



-0.072i far 



and 



leading to the single-particle ANC for the virtual np state 
b = 0.268 fai-V 2 . 



C. The resonance states of 15 F(l/2+, 5/2+) 

Several articles were published recently 0, [2f| [2(| [53, 
l28j testing the predictive power of the current theoreti- 
cal approaches to describe the lowest broad levels in 15 F. 
The final goal of these analyses was a comparison of the 
predictions with the available experimental data on the 
15 F levels. Determination of a broad resonance param- 
eters is a well known unsolved problem in physics. The 
resonance energy and width for a broad resonance are 
not defined uniquely and there are many prescriptions, 
which have been used in literature [29j | . The definitions 
depend not only on the model used, say potential, R- 
matrix, microscopic, but even within a given model the 
pres crip tions for the resonan ce p arameters can be differ- 
ent [29[[3(|. For example, in [3lf four different definitions 
were used. In [29| two more definitions were added in the 
R- matrix approach. That is why we believe that, when 
any compilation includes the broad resonance parame- 
ters, the reference should be done to the prescriptions 
used to determine these parameters. The reason for this 
ambiguity is that for broad resonances in the physical re- 
gion the nonresonant contribution becomes comparable 
with the resonant one. In this case the determined reso- 
nance energy and width depend on how much of the back- 
ground is included into the resonant part. The only way 
to determine correctly the resonance energy and width 
is to single out the resonance pole explicitly in the func- 
tion fitting the experimental data. It is realized in the 
S-matrix pole method. 

Here we address two approaches based on the definition 
of the resonance energy Er = Eq — iT/2 as the energy 
at which the S-matrix has a pole on the second energy 
sheet (low half of the momentum plane) : the potential 
approach based on the solution of the radial Schrodinger 
equation and the analytical expression for the S-matrix. 
The first one gives the most accurate definition of the 
resonance energy and width within the potential model, 
while the second one even more general because it based 
only on the analyticity and the symmetry of the S-matrix 

s 

We remind that a resonance corresponds to the pole 
of the S-matrix at k R — fco ~ i kj and is located in the 
fourth quadrant of the momentum complex plane. Cor- 
respondingly the resonance energy is 



1.2 

2fi 



En 



where 



En 



k 2 



2fi 



(16) 



(17) 



r = 



2 fco fcj 



(18) 



For broad resonances fc/ becomes comparable with fco or 
even larger (fco < fcj). If fcj > fco, i.e. the resonant pole 
in the complex momentum plane, due to large fcj, is far 
from the real energy axis and the energy of the broad res- 
onance, Eq < 0, is located in the third quadrant on the 
second energy sheet and we call it the subthreshold broad 
resonance [5a]. Due to large fcj (or resonance width T), 
the impact of the resonant pole on the cross section or 
scattering phase shift is weakened and the non-resonant 
amplitude or phase shift (non-resonant background) be- 
comes important. The general expression for the elastic 
scattering S'-matrix based on its analyticity and symme- 
try in a vicinity of a single resonance can be written as 



S(fc) 



iS(k) = 2iS p (k) (k - k* R ){k + kg) 

(k - k R )(k + k* R ) 

_ 2i{S p (k)+S R (k)+S a (k)) 



(19) 



where <5 p (fc) is the non-resonant scattering phase shift, 



<5r(&) = — arctan 



fcj 



k — fco 
k — fco 



= — [ arctan 

L 2 ki ' 

is the resonant scattering phase shift [571 ] . and 



Sa(k) 



arctan ■ 



k. 



k + k. 



o 



(20) 
(21) 

(22) 



For narrow resonances, fcj << fco, the phase shift 
|<5 Q (fc)| << 1 can be neglected. In this case, the stan- 
dard method, which we call the phase shift method (or 
"5 = 7r/2" rule), entails the resonance energy Eq the 
value at which the scattering phase S(k) passes through 
7r/2. The resonant width is evaluated from the formula 
r = 2/(dS/dE) at E — Eq or as the energy interval cor- 
responding to change of S from 7r/4 to 3tt/4. However, 
for broad resonances 6 a (k) cannot be neglected and the 
total non-resonant scattering phase shift S p (k) + S a (k) 
becomes dependent on the resonant parameters. This 
non-resonant scattering phase shift may be a large neg- 
ative so that the total phase shift S(k) cannot reach ir/2 
at fc = fco making the n/2 method non-applicable. When 
calculating the elastic cross section or scattering phase 
shift in the presence of the broad resonance, due to the 
importance of the non-resonant phase shift, the cross sec- 
tion depends not only on the resonance parameters Eq 
and F but also on the potential adopted. 

Here as a test case we select resonances representing 
the ground state l/2 + and the first exited state 5/2+ in 
15 F. The latest very detailed analysis of the angular dis- 
tributions for the 14 C(d,p) 15 C reaction [33| showes 
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that the spectroscopic factors for the ground l/2 + and 
the first excited state 5/2+ are close to the single particle 
ones (0.99 and 0.90 correspondingly [33|]. One expects the 
similar numbers for the mirror states in 15 F. Therefore, 
the potential approach is appropriate to describe these 
states. In (34[ the Woods-Saxon potential parameters to 
describe the excitation energies of the mirror levels in 15 C 
and 15 F as well as the experimental data on resonance 
14 +p were found. The authors [34| presented the final 
data on the resonance parameters for the first two levels 
in 15 F using the calculations of the wave function inside 
the nucleus, at the radius of 1 fm. The energy at which 
the absolute value of the wave function reaches its maxi- 
mum was identified as the resonance energy. We call this 
the l^maa;! method. In [3J| the width of the resonance 
was defined by the energy interval over which the ampli- 
tude falls by V2 relative to the maximum of the | ^ m ax \ ■ 
For comparison, in [34j some results were presented using 
also the 7r/2 method. 

We apply the potential 5-matrix pole method by solv- 
ing the Schrodinger equation with the Woods-Saxon po- 
tential given in [34[ for both 15 F resonance states with the 
J n = 1/2+ and 5/2 + . We search for the complex energy 
at which the coefficient c[ ' = (sec Eq. ([7])) similar to 
the search for the bound or the virtual state. Wc note 
that in the standard approach the scattering wave func- 
tion is calculated at real energies, where the non-resonant 
contribution is significant for broad resonances, while the 
Gamow wave function is calculated at the complex en- 
ergy corresponding to the resonant pole of the ^-matrix 
located on the second Riemann energy sheet. As a first 
approximation, to determine the complex resonance en- 
ergy E R = Eq — ir 1 /2 we use the phase shift method 
(or the Ismail method when the 5 = ir/2 method is non- 
applicable). After that, we solve the Schrodinger equa- 
tion near the complex energy E R = Eq — iT 1 /2. The 
final result of this search is the complex energy En, at 
which the coefficient of the incoming wave vanishes. We 
also applied the S-matrix pole search using the analyti- 
cal representation (flH|) for the S-matrix (see explanation 
below) . 

Our results for the energies and widths of the reso- 
nance states are given in Table U compared with the pre- 
vious results obtained using the S = ir/2 and |^ moa ;| 
methods [Hj|. The position E and the width Y of the 
broad resonance depend on the calculation method: the 
5-matrix pole method gives the values of the resonance 
energy and width smaller and more accurate than the 
5 = it/2 and l^ma^l methods. It is worth noting that 
the corrected value of 1.227 MeV for the resonance en- 
ergy of the ground state of 15 F is very close to the lower 
limit given by Fortune obtained using the isobaric 
multiplet mass equation. Besides, in the most recent ex- 
perimental work on 15 F [35[ it was indicated that the 
ground state energy of 15 F can be even lower. 

Figs. [31 and Q] show the real and imaginary parts of the 
normalized Gamow wave function for the l/2 + and 5/2 + 
resonance states in 15 F. The solution of the Schrodinger 



TABLE I: Energy and width of the resonances for the F 
states with J p — l/2 + (the ground state) and 5/2 + (the first 
exited state) calculated by the use of three different methods 
(see the text) 



J p Eq (MeV) T (MeV) Method 



1.450 1.091 
1.290±g; | 0.7 

1.198 0.530 

1/2+ 1.194 0.531 

1.400 0.700 

1.315 0.679 

1.274 0.510 



2.805 0.304 
2.795 ±0.045 0.298 ± 0.06" 
5/2+ 2.780 0.293 
2.777 0.286 

2.762 0.297 



5 = tt/2 

| ^ max 

Pole of S-matrix (potential) 
Pole of S-matrix using Eq. ([19} 
_R-matrix (from 
the scattering phase shift) 
J?-matrix (from the excitation 
function, ro = 4.5 fm) 
J?-matrix (from the excitation 
function, ro = 6.0 fm) 
5 = 7r/2 

^ max 

Pole of S-matrix 
i?-matrix (from the excitation 
function, ro = 4.5 fm) 
_R-matrix (from the excitation 
function, ro = 6.0 fm) 



a It was misprint T = 0.325 MeV for the state | + in Ref. [3. 



equation coincides with the outgoing wave outside the 
potential area. We conclude that the probability of find- 
ing the proton inside the potential area is relatively high. 
The advantage of our method is that we directly find the 
complex energy of the resonant state (energy and width) 
by the same procedure as for the bound state. 

An important test of our method is comparison of the 
single-particle ANC determined as an amplitude of the 
tail of the normalized Gamow function with the ANC 
determined from the residue of the scattering amplitude 
at the pole corresponding to the resonance. For the nor- 
malization of the Gamow wave function we use here the 
method suggested by Zeldovich the numerical appli- 
cation of which is difficult for a broad resonance. How- 
ever, the same relationship between the squared single- 
particle ANC and the residue can be used for both the 
bound and resonance states. One can use Eq. (fl~5"|) 
to find the single-particle ANC of the resonance wave 
function. The results of the calculated residues are 
(-0.038 + i 0.008) fm" 1 and (0.015 - i 0.009) fm" 1 for 
the 1/2+ and 5/2+ states, respectively. From Eq. fp~5|) 
we get the single-particle ANCs (-0.123+3 0.153) fm" 1/2 
and (0.115 + i 0.067) fm" 1 / 2 for the same states, corre- 
spondingly. We obtained the same single-particle ANCs 
directly from the tail of the normalized Gamow wave 
functions validating Eq. (|T5l) . 
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1. Model-independent determination of the energy and 
width of the broad resonance | + in 15 F 

The limitations of the potential model and the exis- 
tence of the phase-equivalent potentials calls for a cross 
check of the energy and width for the broad resonance 
determined from the potential approach. We demon- 
strate how to determine these resonance parameters us- 
ing the model-independent representation of the elastic 
scattering S'-matrix given by Eq. (|19p. Since the ex- 
perimental 2si/2 phase shift for 10 C + p scattering in 
the resonance energy region is not available, we gener- 
ate the "quasi-experimental" 2si/ 2 phase shift using the 
Woods-Saxon potential from [3J], which reproduces the 
14 O + p resonance scattering. Its geometry is ro = 1.17 
fm, a = 0.735 fm, r c = 1.21 fm and the depth V = 53.52 
MeV. The phase shift is shown in Fig [5] Using the S- 




5 10 15 20 25 

r(fin) 




r (fm) 



FIG. 3: Real (a) and imaginary (b) parts of the wave function 
of the l/2 + resonance state in 15 F. The solid line is the 
solution of the Schrdinger equation, the dashed line is the 
outgoing Coulomb function (the Whittaker function). 



matrix pole method from the solution of the Schrodinger 
equation we find the resonance energy for this potential 
E = 1.198 MeV and the resonance width T = 0.530 
MeV. Now we demonstrate that using Eq. (fll)|) we can 
fit the "quasi-experimental" phase shift and determine 
the resonance energy and width. The potential phase 
shift in Eq. (I19|) is approximated by the polynomial 

3 

5 p (k) — ^2 b n (k — k s ) . So, we have 6 fitting param- 

eters including 4 coefficients b n , Eq and T. The final 
result does not depend on the choice of the center of the 
Taylor expansion k s and practically not sensitive to the 
starting values of Eq and V. We take here the starting 
values k s = 0.25 far 1 , E = 1.45 MeV and T = 1.276 
MeV obtained from the 8 = n/2 method, Table [0 The 
fit to the " qausiexperimental" phase shift gives the final 
resonance energy E = 1.194 MeV and T = 0.531 MeV 
what is in a perfect agreement with the results obtained 
using the potential S'-matrix pole method. For the start- 
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FIG. 4: The same as in Fig. [3] but for the 5/2 + resonance 
state in 15 F. 
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ing search values Eq = 1.6 Mev and T = 1.276 MeV we 
get the fitted energy E = 1.198 MeV and T = 0.532 
MeV. Thus Eq. (fP9|) allows one to obtain the energy 
and width of the broad resonance using, for example, 
as input parameters the resonance and width obtained 
by the S — it/2, \*ff ma x\- The model-independent result 
obtained from Eq. (fTO)) gives very close values to the 
potential S-matrix pole. Assigning a 10% uncertainty 
to the "quasi-experimental" phase shift results in a sim- 
ilar uncertainty in the determined "quasi-experimental" 
phase shift resonance energy and width. 



2. Comparison with R-matrix approach 

The resonant S'-matrix obtain from the i?-matrix con- 
tains the nonresonant contribution through the energy 
dependence of the level shift and resonance width. The 
extrapolation of this functions to the complex energy 
plane make them complex, i.e. they lose it physical mean- 
ing. Thus the i?-matrix approach is not designed for ex- 
trapolation to the resonant pole. 

Here we apply the i?-matrix approach to determine 
the energy and the width of the resonance with the S'- 
matrix pole method. For an isolated resonance in the 
single-level, single-channel i?-matrix approach with the 
zero boundary condition the Coulomb-modified nuclear 
scattering S-matrix is 



S 



„2iS h 



Eo-E + i 



r(B) 



E - E -i 



r(E) 



(23) 



where Sh s is the hard-sphere scattering phase shift. To 
obtain this equation the linear energy dependence of the 



T3 




FIG. 5: The 14 0+p 2s 1/2 scatterin; 
the Woods-Saxon potential from [3f 
experimental phase shift" . 



; phase shift generated by 
] and used as the "quasi- 



level shift function A(E) is taken into account [3y]. Here, 
E is the real part of the resonance energy. In the R- 
matrix Eq is determined as E\ + A(Eq) = Eq, E\ is 
the i?-matrix level energy, T(E) = 2 ~/f Pj 2 (E , r ) is the 
observable resonance width in the i?-matrix approach de- 
pending on the energy and the channel radius tq, 7; is 
the observable reduced width amplitude, Pi(E,ro) is the 
penetrability factor in the Z-th partial wave. The res- 
onance width in the i?-matrix approach, in contrast to 
the Breit-Wigner equation, depends on the energy. This 
dependence reflects the fact that the S-matrix in the R- 
matrix is richer than the Breit-Wigner equation: it in- 
cludes also the non-resonant background, which is con- 
tributed by the hard-sphere phase shift and the energy 
dependence of the level shift function and the resonance 
width. For narrow resonance (T(Eq) << Eq) the pole in 
Eq. ([23]) E R w Eq - iT(E )/2. For a broad resonance 
this resonance energy is not a pole of the S-matrix. The 
equation for the resonant pole in this case is given by 
E R = Eq -iT(E R )/2. At complex E R T{E R ) becomes 
complex and loses its meaning of the width. For a broad 
resonance in the i?-matrix method the resonance energy 
is defined as E R = Eq — i T(Eq), which is not a pole of Eq. 
(|23| . Hence, for broad resonances the difference between 
the resonance energy from the S-matrix pole method and 
the i?-matrix method is expected. 

To compare the results for the i?-matrix and S-matrix 
pole methods for the S1/2 4 " resonance we use the phase 
shift generated by the Woods-Saxon potential from (34| 
as the "quasi-experimental" one and determine the res- 
onance energy and width by fitting the i?-matrix phase 
shift to the " quasi-experimental" . The results are shown 
in Table |U The i?-matrix resonance energy and width 
found at tq — 5.0 fm are higher than the S-matrix pole 
ones and close to the l^ma^l result. Both i?-matrix and 
I'I'mazI methods determine the resonance energy from the 
data at real energies where for broad resonances the con- 
tribution of the background becomes important. The 
S-matrix pole method determines the resonance energy 
and width by extrapolating the data to the pole in the 
complex energy (momentum) plane. In the vicinity of the 
pole the resonant contribution becomes dominant com- 
pared to the background and determination of the reso- 
nance parameters is more accurate than in the physical 
region. 

We made additional comparison of the i?-matrix ap- 
proach by fitting the measured in L 34| the excitation func- 
tion of the 14 + p scattering at 180°. Both resonances 
• s i/2 + an( i <^5/2 + coherently contribute to the excitation 
function. The resonances can be separated only after in- 
tegration over the scattering angle. The selection of 180° 
scattering angles minimizes the Coulomb scattering ef- 
fects and enhances the d£, 2 resonance contribution. The 
two-level i?-matrix fitting to the excitation function gives 
the observable resonance energy and width presented in 
Table|I]for two channel radii tq — 4.5 and 6 fm. The reso- 
nance energy is determined as the peak of the \S(k) — 1| 2 , 
and the width as the FWHM of this function. We note 
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that this prescription differs from two prescriptions used 
in [3(3|. For narrow d§/ 2 + all methods gives very close 
results, but it is not the case for the broad resonance 
Si/2 + '. The i?-matrix results are between the l^maxj and 
the S'-matrix pole. Since the S-matrix pole method based 
on Eq. (fT9|) correctly takes into account the resonance 
contribution as a pole in the complex energy (momen- 
tum) plane and analytically continue it to the physical re- 
gion, it allows one to separate correctly the non-resonant 
(background) contribution from the resonance one and, 
hence, provides the most accurate determination of the 
resonance energy and width. 



D. The lowest levels in the mirror nuclei 11 Be and 



The light neutron rich nucleus 11 Be is probably the 
most discussed nucleus. The interest to n Be is related 
to the well known inversion of the shell model levels in 
this nucleus. It has the following low- lying states: ^ 

(ground state), and the excited states \ at E x = 0.320 

MeV and f + at E x = 1.778 MeV [37]. The hrst two 
are the bound states while the third one is a resonance. 
As it was mentioned in Q "the lowering of the si/ 2 or- 
bital with respect to the 0d 5 / 2 orbital is expected for a 
simple potential well". The pi/ 2 state belonging to the 
K = 1/2 band has a pretty stable dominantly [421] spa- 
tial symmetry configuration since the next l/2 _ state is 
10 MeV away [Hj]. In [3^] was shown that that the low- 
est p\/ 2 state obtained in a central potential with the 
spin-orbital interaction strongly overlaps with the state 
projected from a Slater determinant of the lowest orbits 
in the Nilsson's model with the same spin-orbital inter- 
action as the shell model and deformation. 

In this work to test our method we apply it for cal- 
culation of the three lowest states S\/ 2l p\/ 2 and e? 5 / 2 in 
n Be and n N nuclei belonging to the multiplet T = 3/2. 
We also estimate the spectroscopic factors for s 1 / 2 and 
d 5 / 2 states using the potential approach leaving aside Pi/ 2 
state, which is not a single-particle [![. 

Different reactions, including the 10 Be(d, p) 11 Be reac- 
tion with the radioactive 10 Be target, were used to obtain 
the spectroscopic factors for the lowest states in 11 Be. As 
a standard procedure, the single-particle neutron wave 
functions in 11 Be are used as the input in the DWBA code 
to get the neutron spectroscopic factors. The obtained 
spectroscopic factors are in the intervals (0.5-0.96) [ibl ] 
and (0.7-0.8) [4l|. But they are model-dependent because 
they depend on the Woods-Saxon potential adopted for 
the neutron bound state in 11 Be, optical potential in the 
initial and final channels of the (d, p) reaction and accu- 
racy of the DWBA to analyze for the deformed 11 Be [42] ]. 
A priori the transfer reactions involving deformed nuclei 
require the codes, which take into account the multi-step 
transfer mechanisms, for example, the coupled channels 
Born approximation available in FRESCO. That is why 



it is difficult to say from the DWBA analysis to what 
extent the three lowest neutron states are single-particle. 

The nucleus n N is the mirror of 11 Be, and it should 
have a similar level scheme. All states n N are unstable 
to proton decay. Therefore, their decay widths directly 
related to their single particle nature. Since the discovery 
of the ground state in n N in 1996 [43|], the lowest levels 
in n N were studied in many works (see [44j and refer- 
ences therein). In this section, we apply the S'-matrix 
pole method to study the broad levels in n N. Simul- 
taneously we attempt to find restrictions on the single 
particle potentials related to the widths and excitation 
energies of the mirror states in 11 Be and n N. 

To determine the single-particles levels in n Be and 
n N, we use the Woods-Saxon plus Coulomb potential 
similar to the ones used in (2f| |4(| E^, |4(| . The parame- 
ters of the potential are fitted to reproduce the energies 
of the low-lying levels in 11 Be. Then we use this nuclear 
potential plus the Coulomb potential to find the mirror 
levels in n N. We apply the pure single-particle approach 
as in [H, SI]. 

The different sets of the potential parameters, which 
were used to fit the lowest levels in 11 Be, are presented in 
Table [TTJ As a starting point, the standard geometrical 
parameters ro = 1.25 fm, a = 0.65 fm of the Woods- 
Saxon potential are used. Then, we vary the depth of 
the central potential Vq to fit the binding energy of the 
ground state Si/ 2 of 11 Be (well-depth procedure). After 
that, we vary the radius ro and the diffuseness parame- 
ter a to fit the binding energy in 11 Be at the fixed depth 
Vq = 57.057 MeV found from the fitting at standard geo- 
metrical parameters. We use the same procedure for the 
P1/2 and d 5 / 2 states. As is seen in Table HTl the adopted 
potential well is shallower for the pi/ 2 state than for the 
ground state, which reflects the inversion of the s and p 
levels. To reproduce the well-known energy of | + reso- 
nance (Er = 1.275 MeV), we use the set of the potential 
parameters determined for the ground state i + of 11 Be 
with addition of the spin-orbital potential. The fact that 
| + state has particle width provides for an additional cri- 
terion for the selection of the potential. As can be seen 
in Table [Til the calculated single-particle widths for this 
state are larger than the experimental values of 100 ± 20 
keV [S3] and 104 ±21 keV (47[. Taking into account that 

r exp — S T S p (^4) 

we can estimate the spectroscopic factor S for this state. 
Here E sp and T sp stand for Eq and T, correspondingly. 
The experimental and theoretical spectroscopic factors 
are in the range 0.45-0.8 [4l|. The spectroscopic values 
in the interval 0.45-0.61 are obtained by comparing the 
data in Table [TTJ and the experimental ones. Taking into 
account the experimental uncertainties of 20%, the high- 
est value of the spectroscopic factor can be ~ 0.73. (To 
decrease the calculated single-particle width one has to 
use a sharper potential (smaller diffuseness) , which seems 
in contradiction with current experimental data and the 
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TABLE II: Energies and widths calculated for low-lying levels 
of n Be by S-matrix pole method. 
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bound 


5 + 
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1.20 
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57.057 


7.131 


1.275 


0.221 




1.22 


0.713 


57.057 


6.222 


1.275 


0.208 




1.25 


0.650 


57.057 


4.743 


1.275 


0.189 




1.27 


0.607 


57.057 


3.671 


1.275 


0.176 




1.29 


0.562 


57.057 
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0.164 



theoretical predictions (see [34| and references therein)). 
Smaller experimental uncertainties in the width of the 
| + state result in stronger restrictions in the potential 
parameters. 

The S-matrix pole calculations for the three states, 
which are all resonances, for the mirror n N nucleus are 
made using the potential parameters for the 11 Be nu- 
cleus by adding the Coulomb potential of the uniformly- 
charged sphere of the radius parameter rc (Eq.Q). The 
results are shown in Tables IIIII and IIVI for two values of 
the radius of the uniformly-charged sphere. 

It is worth noting that in the case of the relatively 
sharp | + resonance, the differences between calculations 
of the resonance energy and the width using the .S-matrix 
pole and the phase shift are relatively moderate, 140 keV 
and 130 keV, correspondingly. However, these differences 
become significantly larger for the broad resonance i in 
11 N, up to ~ 300 keV for the energy and ~ 500 keV for 
the width. We note that for the same set of the potential 
parameters the S-matrix pole method gives energy and 
width smaller than those obtained by the phase shift. As 
for the 2si/2 state in n N, the phase shift never passes 
through 7r/2 in agreement with the earlier observation by 
Barker [29 ] . 

The S-matrix pole method reveals the resonance pole 
for the state i + in n N. To check that we have found the 
pole correctly, we match the logarithmic derivatives of 
the solution of the Schrodinger equation and the Gamow 
function in the asymptotic region. We also check the ra- 
tio of the solution to the Gamow function, which must be 
constant at the asymptotic region. Let us consider the 



TABLE III: Energies and widths calculated for low-lying lev- 
els of X1 N by S-matrix pole method. The Coulomb radius 
ro = 1-1 fm. 



r 


ro 


a 


V 


Vis 


Esp 






(fm) 


(fm) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


i + 
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1.20 
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57.057 
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0.832 




1.22 


0.713 


57.057 
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0.869 
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0.650 
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0.956 
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37.505 


6.0 
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0.985 




1.27 


0.545 


37.505 


6.0 


2.274 


1.014 




1.28 


0.451 


37.505 


6.0 


2.415 


1.035 


5 + 
2 


1.20 


0.753 


57.057 


7.131 


3.653 


0.946 




1.22 


0.713 


57.057 


6.222 


3.699 


0.913 




1.25 


0.650 


57.057 


4.743 


3.772 


0.865 




1.27 


0.607 


57.057 


3.671 


3.823 


0.834 




1.29 


0.562 


57.057 


2.520 


3.877 


0.804 



| level in n N. By averaging the experimental data from 
Refs. [H, EH, El, m , we obtain 3 - 72 ± 0.050 MeV for the 
resonance energy for this level and 0.55 jlg'^ 5 MeV for the 
width. We can conclude from Tables IIIII and IIVI that the 
S-matrix pole method gives for the width ~ 0.85 MeV 
resulting in the spectroscopic factor S — 0.65. The stan- 
dard geometrical parameters vq = 1.25 fm and a = 0.65 
fm provide for a good agreement with the average ex- 
perimental energy for this level. Using these parameters 
(Table ILTT) . one can find the spectroscopic factor of 0.53 
for the mirror state in 11 Be. Assuming that the spectro- 
scopic factors should be the same for the mirror states, 
one can conclude that the average value of the spectro- 
scopic factor S = 0.59 is a characteristics of the single- 
particle structure for the | + state in the mirror n N and 
11 Be nuclei. 

The experimental data for the broad \ resonance 
state in n N need careful consideration because the re- 
sults reported in [H, |H, |49|, [5(| are different due to dif- 
ferent definitions of "energy" and "width" in these works. 
We nevertheless conclude that the resonance energy of 
this state is ~ 2.2 MeV. As in [H] we use Eq. to get 
the width of the \ state in the potential approach. The 
spectroscopic factor 5=0.66 results in the width T = 0.65 
MeV for this state. This spectroscopic factor coincides 
with the shell model prediction for the analog state of 
11 Be and our result is in a good agreement with the one 
obtained in [381 ] . 

All available experimental data [H 0, H3] give higher 
resonance energies of the ^ + state in n N than our calcu- 
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TABLE IV: The same as in Table |ffl] but for the Coulomb 
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lations (see Table IIIII and IIVI) . (We exclude most of the 
mass-transfer data from the consideration because of the 
very low population of the i state in 11 N in these reac- 
tions.) The n N ground state resonance energies (the rel- 
ative 10 C + p energy) are grouped around 1.3 MeV from 
the data 0, |H(J. The most recent study [5(| resulted 
in the value of 1.54 MeV for the resonance energy; the 
experimental widths for the resonance are in the range 
from 0.83 MeV [5(j to 1.4 MeV 0. These experimen- 
tal values were extracted using different approaches. In 
[43L I44I ] the behavior of the single-particle wave function 
inside the 11 N nucleus is used to determine the resonance 
energy (identified as the energy, at which the modulus of 
the wave function calculated at 1 fm reaches maximum - 
l^moasl method) and the resonance width. The i?-matrix 
analysis was used in 50] . Both approaches can not elimi- 
nate a contribution from the non-resonant potential scat- 
tering. Leaving a detailed analysis of the experimental 
data for future studies, we make a crude estimation of 
the spectroscopic factor for the 2s state. To this end we 
average data from [43|, [44[ for the i + resonance getting 
the resonance energ y 1. 30 MeV and resonance width 1.20 
MeV. The analysis |43l . l44j was based on a potential ap- 
proach, i.e. the Woods-Saxon potential was found, which 
allows fitting the excitation functions and angular distri- 
butions for the elastic proton resonance scattering. Using 
the potential parameters from [43j, we apply here the S- 
matrix pole method, rather than the |^ max | method, to 
determine the resonance energy and width. We obtain 

1.102 MeV for the resonance energy and 840 keV for the 

1 + 

resonance width, i.e. the resonance energy of ^ state in 



n N decreases by 200 keV compared to the one adopted 
previously! We now adopt 1.102 MeV as a new "experi- 
mental" resonance energy of the i state in n N. 

We note that the potential found in [43[ does not re- 

1 4- 

produce the experimental binding energy of the 5 state 

and resonance energy of | + state in n Be. Meanwhile 
the potentials given in Tables IIIII and IIVI fit 2s and Id 
states both in n Be and n N. Then, we assume that the 
potential with the standard geometry ro = 1.25 fm and 
a = 0.65 fm in Table ITVl is the "right" one. Note that 
the resonance energy and width obtained for this poten- 
tial are very close to the average resonance energy and 
width shown in Table IIVI for 5 different potentials. This 
potential gives 1.062 MeV resonance energy, which is a 
pure single-particle energy. We observe that this energy 
is ~ 0.04 MeV less than the "experimental" value of 1.102 
MeV obtained for the potential adopted in [43j. This 0.04 
MeV can be attributed to the non-single-particle admix- 
ture to the structure of the ^ + state n N. 

An estimation of the spectroscopic factor can be ob- 
tained from the consideration of the width of the state. 
The ro = 1.25 fm and a = 0.65 fm parameters generate 

0. 95 MeV for the i + state width at the "experimental" 
resonance energy of 1.10 MeV. The ratio of 0.84/ 0.95 
(the "experimental" width/ calculated width) results in 
the spectroscopic factor of 0.88 for the adopted poten- 
tial with the standard geometry. Hence we obtain much 
larger spectroscopic factors for the ground state in n N 
than for the | + excited state. 

As a final remark to this section, it is worth noting 
that the conventional potential approaches, which deter- 
mine the resonance energy and width from the energy de- 
pendence of the phase shift or from the l^maxl method, 
may not give accurate results because of the distortion 
generated by the non-resonant background at physical 
energies. For example, l^moxl may reach a peak in the 
internal region at energy E 7^ Eq. In this sense the S- 
matrix pole is the most accurate method for a given po- 
tential because it determines the resonance energy and 
width by searching the resonant pole at complex energy, 

1. e. separates the resonant contribution from the back- 
ground. We note that the resonance energy determined 
by the S-matrix pole method depends on the adopted po- 
tential. Moreover, the resonance parameters determined 
by the S'-matrix pole method may differ from the ones 
determined from the i?-matrix approach as we have seen 
it for the 15 F case. 



1. Model-independent determination of the energy and 
width of the broad resonance i + in 11 N 

Here we determine the resonance parameters for the 
2s 1/2 resonance in n N = 10 C + p using the model- 
independent representation of the elastic scattering S'- 
matrix given by Eq. (fH)]) . We use the Woods-Saxon 
potential with the standard geometry ro = 1.25 fm, 
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a = 0.65 fm, the Coulomb radial parameter tq = 1.2 
fm and the depth V — 57.06 MeV to generate the 
"quasi-experimental" 2s\/2 phase shift. The resonance 
energy for this potential obtained from the Schrodinger 
equation is Eq = 1.062 MeV and the resonance width 
r = 0.884 MeV (see Table |IY|. To fit this "quasi- 
experimental" phase shift we use the polynomial approx- 

3 

imation S p (k) = J2 b n (k - k s ) n in Eq. ([19]). The S- 

n=0 

matrix pole method based on Eq. (fT9|) gives Eq = 1.057 
MeV and T = 0.880 which agrees extremely well with 
the potential 5-matrix pole. As starting search value in 
fitting the "quasi-experimental" phase we used k s = 0.25 
fm" 1 but the result only depends slightly on the initial 
k s value. 



E. The subthreshold resonances in the 
proton— proton system 

The poles for the antibound state of the singlet 
neutron-neutron or the neutron— proton systems are lo- 
cated on the imaginary axis in the complex momentum 
plane (at energies E nn = —134 keV and E np = —66 keV 
) . In [l8j using the effective-range approach Kok showed 
that in the case of the proton-proton system the pole 
moves to the complex plane, due to the Coulomb bar- 
rier. The ground state pole of the s-wave pp scattering 
amplitude was found in [3 at k pp = (0.0647 - i0.0870) 
fm -1 or E pp = (-140 - z467) keV. The effective-range 
par ameters for the standard expansion were taken from 
}5 1| . Recently, calculations with the same approximation 
were repeated in [52j resulting in k pp = (0.0644 — iO. 0871) 
fm -1 or E pp = (—142 — i465) keV, which is in a good 
agreement with the Kok's result. 

A definition of the renormalized partial amplitude in 
the presence of the Coulomb interaction was given ear- 
lier (see Eq. (3) in [53|). A new corresponding formula 
was derived in [52j for the renormalized vertex constant 
G ren for the virtual decay of a nucleus into two charged 
particles in the effective-range theory. It was applied 
to the pp and pd systems using the standard effective- 
range expansion and the effective-range function with a 
pole, respectively. The value of G^ e „ is real quantity 
for the bound state because the energy is real. In the 
case of the resonance, the energy is complex so G^ en be- 
comes complex. For the pp ground state, the value 
Gren = (0-060 + z0.051) fm was obtained in [52| with 
the effective-range parameters taken from (5lJ ■ The only 
condition which validates these results is the convergence 
of the effective-range expansion near the pole considered. 
It was shown in [18| that the results change only slightly 
when the parameters of form are neglected (P=Q—0). 
The convergence is ensured in the case of the pp sub- 
threshold resonance pole. 

Nevertheless, the effective-range approximation has 
some drawbacks. It gives the partial scattering ampli- 
tude in an analytical form as a ratio of two polynomials. 



As a result, all the amplitude singularities are poles in 
the complex momentum plane. The number n of the 
poles is obviously defined by the maximal degree used in 
the effective range expansion up to fc™, which gives the 
degree of Kok's equation for the position of the pole and 
correspondingly the number of its solutions without the 
Coulomb force. For example, a logarithmic dynamical 
cut of the amplitude in the case of the two-body model 
with the Yukawa potential can not be reproduced in this 
approximation. But it is imitated by a pole located on 
the positive imaginary axis, which is not a bound state 
pole. In the case of the pp system, the situation is simple 
because there is no bound state, so any pole on the pos- 
itive energy axis is an unphysical one. Moreover, the re- 
gion of the validity of the effective-range approximation is 
limited by the condition \k\ < \k ma x\, where the effective- 
range expansion converges. In the potential model with 
the asymptotic V(r) — ► const ■ r v ■ exp (— r/R), the value 
Ifcmazl = 1/(2J?) is the beginning of the dynamical cut 
on the imaginary axes in the complex k plane. In the 
case of charged particles, the number of roots is infinite 
(see [HI). In particular, as noted in [54| . the sequence 
of poles located near the negative imaginary axis can be 
mistakenly identified as virtual (antibound) state poles 
known for the system without the Coulomb interaction. 

Finding the pole by solving the Schrodinger equation 
is the most reliable way to confirm that the pole found by 
a solution of Kok's equation [18] is not a false one. This 
was done in the present paper for the pp system with 
the Yukawa potential. Its parameters are taken from (5l| 
for the singlet np system. We find that k^ p u = (0.064 — 

i0.082)fm _1 (or E pp = (-106.7 - i 435.5) keV ) for the 
pp ground state. After that, we have slightly changed 
the geometric parameter to describe the experimental pp 
scattering length and effective range and the resolved 
Schrodinger equation gives E pp = (—138.16— i 463.14) 
keV. This result almost coincides with Kok's results. The 
pole for nn system is E nn = — 92 keV. The resonance 
wave function contains the outgoing wave in the asymp- 
totic region while the ingoing wave is absent. For the 
normalization of the wave function in this case we can not 
use the Zeldovich procedure because Re(k pp ) < Im(k pp ), 
however, we can find residue at the pole. The residue in 
the pole corresponding to the subthreshold resonance is 
Ap P = (-0.021 + i0.057) fm" 1 . 



IV. CONCLUSION 

In the present work we apply the S'-matrix pole method 
to determine the energies of the bound, the virtual states 
and resonances. This method is based on a numerical so- 
lution of the Schrodinger equation. Usually this method 
is applied to the bound states, but here it is extended 
to the resonance and virtual states despite the fact that 
the corresponding solutions increase exponentially when 
r — > oo. The method turns out to be especially useful for 
broad resonances including subthreshold ones. 
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There can be a few poles in the complex plane when 
applying the effective range theory and the number of 
the poles is defined by the maximal power k n used in the 
effective range expansion. An additional investigation 
should be done to select a physical pole. In our approach 
one gets no false poles, thus resolving ambiguity problem 
appearing in the effective range approach. 

In the case of the resonances potential models and 
i?-matrix approach are commonly used to analyze the 
experimental data, and resonance parameters are deter- 
mined from the fits. For narrow resonances both ap- 
proaches give accurate results. However, this is not the 
case for broad resonances. In this case, due to the distor- 
tion caused by the non-resonant background at physical 
energies, the resonance energy and the width determined 
from the fitting of the experimental data depend on the 
model and within a given model the prescriptions to de- 
termine the resonance energy and width may be different. 
Usually researchers use different definitions of the reso- 
nance energy and width. Broad resonance parameters 
extracted from the experimental data are model depen- 
dent. For this reason, one should indicate the method 
used to determine them in any subsequent references. 

Here we address two methods for determining the res- 
onance energy and width from the pole of the S-matrix: 
the potential S'-matrix pole method based on the solu- 
tion of the Schrodinger equation and the S'-matrix pole 
method based on the analytical continuation for the S'- 
matrix to the resonant pole. We compare the results 
for the resonance parameters obtained from the differ- 
ent determinations of the resonance energy and width in 
the potential approach, the S-matrix pole methods and 



_R-matrix method. Correct evaluations of the resonance 
parameters are important when comparing the experi- 
mental data, both for the tests of the isobaric multiplct 
mass equation and for detailed structure calculations of 
the exotic nuclei. The potential S-matrix pole method 
provides the most accurate resonance energy and width 
for a given potential. The second S-matrix pole method, 
which uses Eq. 1191 is even more general because it does 
not require any potential model and is based only on the 
analyticity and symmetry of the S-matrix. In contrast to 
other approaches, the pole S-matrix methods allow one 
to correctly separate the resonance pole contribution and 
the nonresonant background. 

Our approach has a potential of being extended to 
treat broad resonance populated in transfer reactions, 
where the half-off-cnergy shell resonant amplitude inter- 
feres with the half-off-energy shell nonresonant ampli- 
tude. At present there is a huge disagreement in the res- 
onance parameters for broad resonances obtained from 
the resonance or direct reactions 1441. 
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